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ABSTRACT 

This  reporC  presents  refinements  of  a previous  momentum- 
integral  method  for  calculating  three-dimensional  turbulent 
boundary  layers  on  ship  hulls.  In  particular  the  following  refine- 
ments are  made:  the  small  crossflow  assumption  is  removed; 
numerical  calculation  of  the  double  model  potential  flow  replaces 
the  slender  body  potential  flow;  a more  general  and  versatile  or- 
thogonal coordinate  system  is  used  in  place  of  the  streamline 
surface  coordinate  system;  and  finally,  an  improved  numerical 
method  is  used  for  solving  the  momentum- integral  boundary-layer 
equations.  It  is  shown  that  the  boundary  layer  calculation  method, 
developed  here,  can  be  used  to  calculate  certain  boundary  layer 
parameters,  such  as  boundary  layer  thickness  or  skin  friction, 
with  fair  accuracy  over  a large  portion  of  hulls  that  maintain  un- 
separated flow.  The  surface  coordinate  system  can  also  be  used  in 
other  methods  for  calculating  the  boundary  layer. 


ADMINISTRATIVE  INFORMATION 

The  work  reported  herein  was  supported  by  the  General  Hydromechanics 
Research  Program  under  Task  Area  SR  023-01-01  and  Work  Unit  1552-070. 


INTRODUCTION 

This  report  presents  some  refinements  of  the  momentum- integral  method 

for  calculating  three-dimensional  turbulent  boundary-layers  as  developed  by 
1* 

von  Kerczek  for  ship  hulls.  These  refinements  include:  (1)  the  removal 
of  the  boundary-layer  small  crossflow  approximation;  (2)  the  incorporation 
of  an  exact  numerical  calculation  of  the  double  model  potential  flow  in- 
stead of  using  the  slender-body  theory  potential  flow;  (3)  the  abandonment 
of  the  streamline  surface  coordinate  system  in  favor  of  a more  general  and 
versatile  orthogonal  surface  coordinate  system;  and  (4)  an  improved 
numerical  method  for  solving  the  momentum- integral  boundary  layer  equations. 

The  use  of  momentum- integral  methods  for  calculating  three-dimensional 
boundary  layers  has  come  under  severe  criticism  recently  by  the  advocates 

of  the  differential  boundary  layer  equations  (see,  for  example,  Landweber 
2 3 4 

and  Patel,  Cebeci  et  al.  and  Spalding  ).  The  main  objection  to  the  inte- 
gral methods  seems  to  center  on  the  unavailability  of  a suitable  crossflow 
velocity  profile  function  that  adequately  approximates  a variety  of 


*A  complete  listing  of  references  is  given  on  page  67. 
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crossflow  velocity  profiles.  It  is  most  often  claimed  by  these  critics 
that  an  accurate  representation  of  the  boundary  layer  crossflow  profile  is 
required  for  the  prediction  of  longitudinal  bilge  vortices.  Large-scale 
longitudinal  bilge  vortices  arise  due  to  a complicated  form  of  three- 
dimensional  separation;  the  vortex  flow  itself  being  the  separated  flow. 
Thus,  one  cannot  expect  to  be  able  to  compute  the  bilge  vortex  flow,  even 
by  the  most  sophisticated  boundary  layer  methods,  whether  integral  or 
differential.  Presently,  boundary  layer  theory  can  be  used  only  to  calcu- 
late the  flow  up  to  separation.  Recently  developed  momentum- Integral 
methods  for  two-dimensional  and  axisymmetric  boundary  layers  (see 
Green  et  al.^)  are  as  accurate  as,  yet  considerably  more  economical  than, 
differential  methods.  There  seems  to  be  no  reason  to  believe  that  similar 
improvements  in  three-dimensional  momentum- integral  methods  cannot  be  found. 

Present  models  of  the  boundary  layer  crossflow  are  primitive  and 
further  experimental  data  and  research  can  be  expected  to  uncover  a simple 
crossflow  velocity  profile  family  that  is  adequate  for  calculation  methods. 
For  this  reason  it  was  thought  desirable  to  make  the  technical  improvements 
mentioned  above  in  the  von  Kerczek  method  so  as  to  accommodate  easily 


detailed  improvements  in  the  crossflow  model  that  may  come  about  later. 

However,  the  modified  surface  coordinate  system  and  the  incorporation  of 

the  exact  double  model  potential  flow  calculations  (in  lieu  of  the  slender 

body  theory  potential  flow  calculation  method)  are  of  independent  value  and 

can  be  used  with  any  other  boundary  layer  calculation  method.  The  surface 

coordinate  system  for  the  boundary  layer  calculations  developed  in  this 

report  has  some  advantages  over  the  coordinate  systems  recommended  by 

others  (see,  for  example,  Ccbect  et  al.  and  Miloh  and  Patel  ) . The  surface 

coordinate  system  used  in  this  report  is  very  similar  to  the  one  used  by 
3 

Cebeci  et  al.  but  it  does  not  have  the  complication  of  being  nonorthogonal . 

£ 

The  present  surface  coordinate  system  is  superior  to  the  Miloh  and  Patel 
coordinate  system  because  it  provides  a better  coordinate  net  coverage  of 
the  hull  surface  for  uniform  spacing  of  the  coordinate  parameters. 

This  report  is  divided  into  six  sections  including  the  introduction. 

The  second  section  describes  the  formulation  of  the  boundary  layer 
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calculation  problem  in  terms  of  the  momentum- integral  entrainment  method. 

The  third  section  describes  the  surface  coordinate  system  and  the  potential 
flow  calculation  method.  The  fourth  section  describes  the  numerical  method 
for  integrating  the  momentum-integral  equations.  The  fifth  section 
describes  and  discusses  some  sample  computational  results  and  the  sixth 
section  gives  some  concluding  remarks  on  further  developments  of  this 
three-dimensional  ship  boundary  layer  calculation  method.  An  appendix  at 
the  end  of  the  report  gives  some  detailed  formulas  that  are  used  in  the 
computational  algorithm. 

THE  BOUNDARY  LAYER  MOMENTUM  INTEGRAL  EQUATIONS 
It  is  assumed  that  the  ship  surface  is  hydraulically  smooth  and  has 
no  abrupt  changes  in  principal  curvature  anywhere  on  the  hull.  Also  assumed 
is  that  the  boundary  layer  thickness  is  small  compared  to  the  principal 
curvatures  everywhere  on  the  hull.  Only  the  turbulent  boundary  layer 
development  is  considered.  The  length  scale  used  is  half  the  length,  L, 
between  the  perpendiculars  of  the  hull.  The  velocity  scale  is  the  steady 
ship  speed  U^.  Henceforth,  all  physical  quantities  that  are  discussed  will 
be  dimensionless  with  respect  to  these  scales.  There  is  an  orthogonal 
surface  coordinate  system  on  the  hull,  which  has  lines  of  constant  8 running 
generally  lengthwise  along  the  ship  and  lines  of  constant  <p  running  nearly 
parallel  to  the  cross-sections  of  the  ship.  This  system  will  be  described 
in  detail  in  Section  3.  The  coordinate  perpendicular  to  the  hull  surface 
can  be  described  in  terms  of  its  arc  length  parameter  X. 

At  an  arbitrary  point  on  the  hull  surface  the  potential  flow  velocity 
vector  U is  given  by 


U 


u*  + ue  Se 


(l) 


where  e^  and  e^  are  unit  tangent  vectors  in  the  direction  of  the  <J>  and  0 
coordinates,  respectively,  and  U ■ ||u||  is  the  magnitude  of  the  velocity 
vector  U.  The  angle  that  the  velocity  vector  U makes  with  the  <p  coordinate 
line  is  denoted  by 


3 


In  terms  of  this  coordinate  system  the  momentum- integral  boundary-layer 

equations  are  a special  case  of  the  equations  given  by  Myring^  and  repro- 

8 

duced  in  the  article  by  Reynolds  and  Cebeci.  These  equations  are  the 
momentum  integrals  for  the  flow  in  the  ({(-direction; 
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coordinates,  respectively;  0^,  0^,  0^,  0„9,  A^,  and  A^  are  momentum  and 
displacement  thicknesses  defined,  respectively,  by  Equations  (4a, 4f);  and 
C.  and  C.  are  wall  skin  friction  coefficient  components  in  the  $ end  6 


directions,  respectively.  The  usual  turbulent  boundary  layer  assumptions 
of  the  neglect  of  turbulent  normal  stresses  and  the  neglect  of  mean  dif- 
fusion in  directions  parallel  to  the  hull  surface  are  incorporated  in 
Equations  (3a, 3b). 

The  momentum  and  displacement  thicknesses  are  defined  by 


where  u and  v are  the  components  of  the  mean  boundary  layer  velocity  in  the 
<f>  and  6 directions,  respectively. 


In  the  case  of  streamline  coordinates,  in  which  the  coordinate  curves 
♦ are  parallel  to  the  inviscid  streamlines  on  the  hull  surface,  U0  = 0 and 
U$  " U*  For  thi«  case,  the  momentum  and  displacement  thicknesses  of 
definitions  (4a-4f)  then  will  be  denoted  by  corresponding  lowet  case  Greek 
letters  and  these  definitions  reduce  to 


(5a) 


(5b) 


(5c) 


(5d) 


(5e) 


(5f  ) 


The  quantity  6 is  some  overall  nominal  boundary  layer  thickness.  The 
relationships  between  the  boundary  layer  thicknesses  defined  by  Equations 
(4a-4f)  and  those  defined  by  Equations  (5a-5f)  will  be  needed  later  and 
can  readily  be  worked  out.  They  are 


6 


2 2 
©11  ■ 0^  cos  a - (6^2+®21^  sin  a cos  a + ®22  sin  a 


2 2 

©12  ■ (®u“®22^  sin  a COS  a + ®12  COS  Q ” ®21  Sin  a 


2 2 

®21  " ^®h“®22^  sin  “ COS  a + ®21  COS  a ” ®12  sin  a 


2 2 
©22  “ sin  ot  + (612+®21^  sin  a cos  a + ®22  cos  a 


A^  - 6^  cos  a - 62  sin  a 


A2  “ 6^  sin  a + ($2  cos  a 


The  inverses  of  Equations  (6a-6f)  that  relate  the  lower  case  Greek  symbols 

6 and  6 to  their  upper  case  counterparts  easily  can  be  derived. 

The  two  momentum  equations,  Equations  (3a, 3b)  are  insufficient  to 

determine  the  eight  unknowns  0^,  ©12>  ©21>  ©22>  A^  A2,  Cf  and  Cf  , hence, 

Q <p 

some  other  integral  equations  and  empirical  information  are  needed.  The 

9 

calculation  method  is  based  on  the  Cumpsty  and  Head  momentum- integral 

8 7 

method  which  utilizes  the  three-dimensional  entrainment  equation 


cy-nij)  - k+  cy-uv 


+ db  <V-UV  - Ke  <V-uvl  ' F 

0 


where  F is  the  three-dimensional  rate  of  entrainment  function.  It  is 
assumed,  in  this  three-dimensional  boundary  layer  calculation  method,  that 


F Is  Che  same  function  evaluated  with  respect  to  the  flow  components  in  the 
streamline  direction  as  the  one  used  for  two-dimensional  flows.  Thus, 
according  to  Green,  et  al. 


F(H)  E 0.025  H - 0.022 


where 


Furthermore,  it  is  assumed  that 


_ *1 
H = 7T- 
911 


6 E 011(G+H) 


where 


G = 2 H 

' (H-l) 


Thus,  the  entrainment  parameters,  the  length  scale  6 and  the  entrainment 

function  F,  are  completely  specified  in  terms  of  integral  thicknesses 

defined  by  Equations  (5a-5f).  A justification  for  the  foregoing  empirical 

q 

expressions  can  be  obtained  from  the  paper  of  Cumpsty  and  Head  and  its 
antecedents. 

The  components  C,  and  C,  of  the  skin  friction  coefficient  can  be  ob- 

talned  from  the  skin  friction  coefficient  in  the  local  inviscid  stream- 
line direction  by  the  formulas 


■ C^(cos  a - sin  a tan  3) 

$ 


(11a) 


Cf  « C,(sin  a + cos  a tan  8) 

e 1 


(lib) 


where  6 is  the  angle  between  the  direction  of  the  wall  friction  vector  and 
the  inviscid  streamline.  The  angle  3 is  precisely  defined  by 


3v 


tan  0 ■ limit  — ■ limit 

**  u ^ lx 


(12) 


where  the  boundary-layer  velocity  components  u and  v are  in  the  direction 
of  the  inviscid  streamline  and  its  normal,  respectively.  By  assuming  that 
the  component  of  the  boundary  layer  flow  in  the  inviscid  streamline  direc- 
tion satisfies  the  two-dimensional  velocity  similarity  laws,  the  skin 
friction  coefficient  can  be  evaluated  using  a two-dimensional  skin 
friction  formula.  The  following  skin  friction  formula  given  by  Head  and 


Patel^  is 


■ exp(aH+b) 


(13) 


where  a - 0.019521  - 0.386768  c + 0.028345  c2  - 0.000701  c3 


b - 0.191511  - 0.83489  c + 0.062588  c2  - 0.001953  c3 


In  R, 
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and  where 


'11 


•m 


ue 


li 


(14) 


is  the  local  streamline  momentum  thickness  Reynolds  number.  Recall  tha*-  U 
is  the  dimensionless  (scaled  by  U^)  magnitude  of  the  local  inviscid  velocity 


at  the  edge  of  the  boundary  layer  and  6^  1b  the  streamline  component  of 
the  momentum  thickness  made  dimensionless  by  the  half  length  L of  the 
ship. 

The  final  empirical  formulas  that  are  used  to  reduce  the  number  of 
unknown  quantities  to  three,  in  order  to  match  the  number  of  differential 
Equations  (3a, 3b)  and  (7),  are  the  formulas  that  relate  the  crossflow 
momentum  and  displacement  thicknesses  012,  021,  ©22,  and  62  to  the  stream- 
line momentum  and  displacement  thicknesses  0^  and  6^.  These  relationships 
constitute  the  critical  approximations  of  three-dimensional  boundary  layer 
theory  in  the  momentum- integral  framework.  It  is  easiest  to  simply  make 
crossflow  and  streamline  flow  boundary-layer  profile  assumptions  and  derive 
the  corresponding  relationships  that  result  from  the  definitions  (5a-5f). 
However  it  is  not  really  necessary  to  proceed  in  this  way.  Basically, 
definitions  (5a-5f)  simply  say  that  if  the  streamwise  velocity  profile  u is 

and  the  crossflow  pro- 

i**l, . . . ,n 

then  the 

j-1, . . . ,m 

boundary  layer  thicknesses  0^  and  6^,  k and  H - 1,2  are  each  functions  of 
the  parameters  [c^],  [ 3 j ] such  as  0^  (o^, . . . ,c»n;  3^,...,^).  One  could 
determine  these  functions  experimentally  and  thereby  completely  bypass  any 
velocity  profile  assumptions. 

In  fact,  such  a scheme  has  already  been  used  for  the  streamwise  flow 
by  using  the  Head  entrainment  method.  From  two-dimensional  and  axisymmetric 
flow  theory,  the  entrainment  method  gives 


described  parametrically  by  parameters  [a^J 
file  v is  described  parametrically  by  the  parameters  [3j] 


«1  " (S1(H,011) 


(15a) 


Cf  - Cf(H,011) 


(15b) 


directly  without  any  profile  assumptions.  In  the  most  highly  developed 
entrainment  method  of  Green  et  al.,^  a third  Independent  parameter  is  added 
to  0^  and  H,  namely  the  entrainment  coefficient  C£,  so  that 
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{1  ■ {1<H'6U-CE> 


<16) 


replaces  Equation  (15a).  If  a velocity  profile  la  needed,  then  the  velocity 
11  12 

profiles  of  Coles  or  Thompson  may  be  used  for  specified  values  of  the 

momentum  thickness  0^  and  shape  factor  H.  Sufficient  experimental  data 

for  three-dimensional  boundary  layers  are  not  yet  available  for  a similar 

program  to  be  conducted.  It  does,  however,  seem  that  the  flow  in  the 

streamwise  direction  is  sufficiently  similar  to  two-dimensional  flow  that 

the  two-dimensional  data  can  be  directly  applied  to  this  component  of  the 

boundary  layer  flow.  However,  insufficient  data  exist  for  the  crossflow  to 

make  more  than  a crude  estimate  of  profile  shapes.  There  has  not  even  been 

a sufficiently  large  collection  and  analysis  of  crossflow  profiles  to  make 

a reasonable  estimate  of  the  proper  parameters  [0.]  that  need  to  be 

j "1  * ...  i m 

used  to  approximate  the  crossflow.  Thus,  as  the  simplest  first  approxi- 
mation of  the  crossflow  profile  shape,  it  is  common  practice  to  assume  that 
(1)  the  crossflow  profile  scales  on  the  same  length  scale  6 as  the  stream- 
wise  profile  and  (2)  the  crossflow  profile  shape  depends  on  only  one  or  two 
independent  parameters,  one  of  which  is  the  shear  stress  angle  0.  The  fact 
that  condition  (12)  must  be  satisfied  at  the  wall  introduces  the  angle  0 
into  the  description  of  the  crossflow  profile  and  also  dictates  (as  a matter 
of  convenience)  the  shape  assumption  in  the  form 


J - (jj)f(X)  tan  & 


(17a) 


In  order  to  satisfy  condition  (12)  at  the  wall  and  the  condition  that 

v « 0 at  A ■ 6 (in  streamline  coordinates),  the  function  f must  satisfy 

13 


f(0)  ■ 1 and  f (6 ) » 0.  The  assumption  by  MagerJ 
proximation  that  gives 


is  a popular  first  ap- 


f(A) 


(i-  I) 


07b) 
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It  should  be  emphasized  that  a profile  assumption,  such  as  Equation  (17b), 

need  not  be  employed  at  all  but,  at  this  time,  lack  of  experimental  data 

forces  the  use  of  such  an  approximation.  More  sophisticated  profile 

1A 

assumptions  than  Equation  (17b)  have  been  made.  For  instance  Okuno 
assumed  that 


-[MD]  (x-t) 


The  parameter  C in  Equation  (18)  is  governed  by  an  additional  crossflow 

1A 

integral  equation.  Okuno  chose  the  crossflow  moment-of-momentum  integral 

equation  for  calculating  the  development  of  the  parameter  C and  obtained 

considerably  better  results  for  crossflow  profile  predictions  than  with  the 

Mager  model  of  Equation  (17b).  However,  there  still  were  certain  areas  on 

the  ship  hull  at  which  the  predicted  crossflow  velocity  profiles  were  in 

serious  disagreement  with  experiment.  The  areas  of  serious  disagreement 

between  the  crossflow  profiles  given  by  Equation  (18)  and  the  experimental 

ones  on  the  Okuno  test  model  (a  Series  60  block  0.70  double  model)  are  very 

close  to  the  stern  on  streamlines  that  turn  upwards  from  the  keel  towards 

1A 

the  load  waterline.  From  Okuno' s experience,  it  seems  likely  that 
further  experiments  and  research  will  eventually  lead  to  a fairly  accurate 
crossflow  velocity  profile  shape  function  f (X ) that  involves  only  one  or  at 
most  two  extra  parameters. 

This  report  is  mainly  concerned  with  setting  the  proper  framework  for 
the  three-dimensional  momentum- integral  boundary  layer  computational  method, 
so  it  will  presently  be  confined  to  the  simplest  of  the  crossflow  models, 
namely  that  of  Mager,  Equation  (17b).  It  is  hoped  that  future  developments 
and  availability  of  sufficient  experimental  three-dimensional  boundary  layer 
data  will  warrant  modifications  of  this  method  to  include  a more  complete 
and  accurate  crossflow  model  along  the  lines  of  the  Okuno  model. 

By  examining  definitions.  Equations  (5b),  (5c),  and  (5f),  it  is  easy 
to  see  that  0j^»  ©2^»  and  62  satisfy  the  relationship 


(19) 


921  ’ *2  * «12 


This  equation  can  be  used  to  eliminate  one  of  the  unknowns  from  the  boundary 
layer  momentum-integral  equations.  The  Mager  crossflow  model.  Equation  (17), 
in  conjunction  with  the  approximation 


u m(H'1)/2 

U “ \&) 


(20) 


and  Equation  (10)  for  the  nominal  thickness  6,  can  be  used  to  express  the 
crossflow  boundary  layer  thicknesses  in  terms  of  the  streamwise  momentum 
thickness  0^,  the  shape  parameters  G and  H,  and  the  crossflow  angle  3 by 
the  equations 


2 6ii^G+H)  tan  B 

021  " H(H+1)  (H+2)  = 911  tan  9 ^1(H) 


(21a) 
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- eu(o+H>  c.  s [i  - 5^  + jjij  ♦ jJj  - 5^]  i eu  t«  6 *2<h)  (22b) 

- 9u(o«)  t«»2  e[f  - + - BJ  + St]  3 °u  t,n2  6 f3<H><2lc> 

s2  ■ - 9u<°**>  t,n  8[ii=r  - iSJ  + Si]  5 011  ta"  6 ?4<">  <21»> 


Thus  Equations  (8)  through  (11),  (13),  (17),  (19),  and  (21),  together  with 
the  transformation  Equations  (6a-6f),  can  be  used  to  reduce  the  total 
number  of  unknown  quantities  to  three.  A convenient  set  of  unknown 
quantities,  that  are  integrated  in  the  (<f>,9)  coordinate  system  by  Equations 
(3)  and  (7)  are  the  streamwise  momentum  thickness  0^,  the  shape  factor  H, 
and  the  tangent  of  the  wall  crossflow  angle  t ■ tan  3.  The  details  of  the 
final  forms  of  Equations  (3)  and  (7)  in  terms  of  these  variables  are  given 
in  the  appendix. 


THE  SHIP  SURFACE  COORDINATE  SYSTEM 

The  hull  surface  coordinate  system  that  is  used  in  the  boundary  layer 
calculation  method  described  in  the  previous  section  stems  directly  from  the 
hull  surface  representation  of  von  Kerczek  and  Tuck.15  This  hull  surface 
representation  utilizes  conformal  mapping  onto  a unit  circle  of  the  cross 
sections  of  the  hull  and  polynomial  interpolation  along  the  length  of  the 
hull  of  the  individual  mapping  coefficients.  Let  s be  the  longitudinal 
coordinate,  x the  lateral  coordinate,  and  y the  vertical  coordinate  of  the 
ship  hull.  Taese  coordinates  are  made  dimensionless  by  the  half  length  L 
so  that  the  bow  and  stern  perpendiculars  lie  at  s - + 1,  respectively.  The 
load  water-line  is  located  at  y - 0 and  the  keel  at  midships  is  located  at 
y - -D. 

The  hull  surface  representation  of  von  Kerczek  and  Tuck15  results  in  a 
parametric  surface  equation  of  the  form 


z 


N M 

x + ly  ■ zz 

n«l  m»l 


. i(3-2n)6  m-i 

A e s 

mn 


(22) 


where  the  matrix  (A^)  of  coefficients  specifies  the  hull  form  and  is 
computed  from  a set  of  defining  hull  offsets  by  an  algorithm  given  by 
von  Kerczek  and  Tuck.15 

The  surface  coordinate  system  used  in  the  boundary  layer  equations 
consists  of  the  9 - constant  lines,  obtained  from  Equation  (22)  and  their 
orthogonal  trajectories,  here  denoted  by  * - constant  lines.  It  is  not 
necessary  to  specify  the  variable  * since  the  arclength  along  the  * - 
constant  lines  will  be  used  directly. 

Equation  (22)  can  be  written  in  real  form  as 


x ■ x(s,0) 


N M 

a'n~1  «s<*-2n)0 

n“ l m-1 


(23a) 
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(23b) 


N M 

y - y(s,0)  - X sm_1  sin(3-2n)6 

n»l  m-1 


and,  in  the  vector  form, 

r « r(s,0)  = x(s,0)i  + y(s,0)}  + sk  (24) 

where  r is  the  vector  from  the  origin  of  the  (x,y,s)  coordinate  system  to  a 
point  on  the  hull  surface,  and  (i»i»k)  are  unit  tangent  vectors  to  the 
(x,y,s)  coordinates,  respectively. 

The  surface  coordinate  lines  0 **  constant,  run  along  the  length  of  the 
hull  surface,  and  the  coordinate  lines  <f>  * constant,  are  nearly  parallel  to 
the  hull  cross-sections.  The  arclength  increments  along  the  <{>  and  0 
coordinates  are  dZ^  and  dZ^,  respectively,  and  are  given  by 


(dr *dr) 


1/2 

©■constant 


• [(H)2+(ii)2«] 


1/2 


ds 


(25a) 


and 


dZ, 


(dr*dr)^2 
- - <Ji=constant 


3x  ds  V 
9s  d0  <f>  / 


+ 


(3y  . 3y  ds 
30  + 3ld0 


1/2 


d0 


(25b) 


where  (ds/d©^)  denotes  the  evaluation  of  the  derivative  ds/d0  along  the 
4>  ■ constant  line. 
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Let  e^,  eg,  and  en  be  unit  tangent  vectors  to  the  <p,  0,  and  the  hull 

surface  normal  coordinates,  respectively.  Then  e„  and  e are  easily 

— b — n 

computed  using  Equation  (24)  by 


3r 


(26) 


e 

-n 


(27) 


where  | | • | | denotes  the  length  of  the  vector.  The  requirement  that  the  <J> 
and  0 coordinates  be  orthogonal  imposes  the  condition  that 


-0 


(28) 


The  increment  of  arc  along  the  curve  <J> 
written  as 


constant. 


can  be 


/ 3r  3r  , \ 
d-|<t>  ~\30  + Ite  df  / d9 


(29) 


and  by  virtue  of  orthogonality 


e,  • dr  , = 0 
-<p  - 4> 


(30) 


The  derivative  (ds/d0)^  can  be  obtained  from  Equation  (30)  by 


ds 

d0 


( 3x 

4.  Iz  h.) 

V 30 

3s 

30  3s/ 

/Sx  \ 

2 

/3v\  2 1 

(30) 

+ 

-2-Z.  I a.  i 

y30 ) + 1 

(31) 
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The  geodesic  curvature  terms  and  are  defined,  respectively,  by 


% " 


d?e 


d*. 


(32a) 


de, 

K„  - er  • — ^ 


0 -0  d* 


0 


(32b) 


The  derivatives  in  Equations  (32a, 32b)  can  be  evaluated  by  converting  them 
to  derivatives  with  respect  to  s and  0 using  Equations  (25a, 25b), 
respectively. 

Miloh  and  Patel^  recommended  the  use  of  the  orthogonal  surface 
coordinate  system  that  consists  of  the  cross-section  curves  and  their 
orthogonal  trajectories  on  the  hull  surface.  In  terms  of  the  surface 
Equation  (23)  the  cross-section  curves  are  given  by  s ■ constant.  Thus, 
if  e,  in  this  case,  is  the  unit  tangent  vector  to  the  cross-section  profile, 
the  orthogonal  trajectories  of  the  cross-sections  can  be  computed  by  inte- 
grating the  differential  equation 


dr  • 0 


(33) 


which,  in  expanded  form,  and  making  use  of  Equation  (29),  reads 


d0 

ds 


9r 

9s 


(s-S) 


(34) 


Examples  of  hull  surface  coordinate  grids  for  the  (<f> .0)  system 
described  earlier  and  the  cross-section  system  of  Equation  (34)  are  given 
in  the  section  on  Computational  Results  and  Discussion. 


1 


17 


The  potential  flow  velocity  on  the  surface  of  the  hull  was  obtained  by 
using  the  computer  program  of  Chang. ^ This  program  solves  the  double- 
model Neumann  problem  by  distributing  a layer  of  doublets  on  the  hull  sur- 
face and  numerically  solving  the  resulting  integral  equation  by  a panel 
method.  The  main  advantage  of  using  a doublet  distribution,  rather  than  a 
source  distribution  as  in  the  Hess-Smith  method,17  Is  that  the  surface 
potential  is  obtained  directly  as  the  solution  of  the  integral  equation. 

The  hull  surface  representation.  Equation  (23),  is  used  to  generate 
the  input  for  the  potential  flow  program  of  Chang.  A uniform  rectangular 
distribution  of  points  in  the  (<t>,6)  plane  determines  a set  of  curved  quadri- 
lateral elements  covering  the  hull  surface.  These  curved  quadrilateral 
elements  are  approximated  by  plane  quadrilateral  elements  and  then  used  as 
input  into  the  Chang  program.  The  results  of  the  potential  flow  calcu- 
lation are  values  of  the  surface  velocity  potential  at  the  geometric  mid- 
points of  the  plane  quadrilateral  elements.  These  values  of  the  surface 
velocity  potential  are  then  assumed  to  be  accurate  values  of  the  exact 
surface  velocity  potential  at  the  points  of  the  hull  surface  that  correspond 
to  the  geometric  centers  of  the  rectangular  elements  in  the  ($,6)  plane. 

The  value  of  the  surface  velocity  at  each  point  is  obtained  by  numerical 
differentiation  of  the  surface  velocity  potential.  The  values  of  the 

potential  are  first  interpolated  along  4>  ■ constant  curves  by  a periodic 
18 

cubic  spline.  This  interpolation  yields  accurate  values  of  the  surface 
velocity  potential  at  arbitrary  locations  on  the  0 = constant  curves.  Then 
the  values  of  the  surface  velocity  potential  are  interpolated  along 
9 » constant  curves  by  another  cubic  spline.  The  surface  velocities  and 
the  derivatives  of  the  surface  velocities  that  are  required  in  Equations 
(3)  and  (7)  are  obtained,  respectively,  by  differentiating  the  cubic 
spline,  evaluating  the  result,  which  gives  the  velocities,  and  then 
refitting  the  velocities  with  cubic  splines  and  differentiating  the  second 
set  of  splines.  This  "spline-on-spline"  procedure  seems  to  be  one  of  the 
best  ways  of  obtaining  two  derivatives  of  a numerically  defined  function. 


NUMERICAL  ANALYSIS 

The  two  momentum  integral  Equations  (3a, 3b),  together  with  the  entrain- 
ment Equation  (7),  can  be  written  completely  in  the  form 


D(W(<t>,0))  If-  + B(W(<$>,6))  |f-  - C(W) 
d\  3£0 


(35) 


where  the  vector  W is  given  by 
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From  Equations  (50)  and  (55)  in  the  appendix,  the  coefficient  matrices  are 
defined  as  follows: 


D(W)  - (D^) 


: $ 

/8 


11 


’21 


11 


0 


9g 


11 


0 


11  9t 


0 


3g 


11 


21 


11  9H 

3g 


11  9t 
9h 


0 


21 


0 


11 


11  9H 
9h 


(37a) 


0 


11  9t 
3g 


0 


11 


12 


B(W) 


(Bij) 


11  9t 

3g 


0 


11  9H 


9g12 
11  9H 


\ 


0 


22 


22  "11  9t 


0 


9g 


22 


11  9H 


3^22  3^22 
h A A ^ 

'22  W11  9t  yil  9H 


(37b) 


/ 


19 


The  variables  in  Equation  (35)  and  auxiliary  formulas  of  Equation  (37)  are 
defined  in  terms  of  the  principal  unknown  momentum  and  displacement  thick- 
nesses 0^,  0^2’  ®21*  ®22’  ^1’  an<*  ^2  *n  t*ie  aPPen<^^x*  *-s  assumed  that 

the  values  of  W are  known  at  a given  station  <J)  = <f>0  for  all  values  of  0 and 
that  the  boundary  layer  is  to  be  computed  between  <f>0  and  a final  station  <p^ 
downstream  of  4>0.  Because  of  symmetry,  it  is  only  necessary  to  solve 
Equation  (35)  between  9 = - n/2  and  0=0.  An  (N+l)  * (Mfl)  grid  with 
spacing  A<p  and  A0,  respectively,  is  superimposed  on  the  region  [ (<j>,0)  j 4>q< 
<}><<}>£, -tt/2<9<0]  . For  simplicity  of  notation5  will  denote  W(<j>Q+iA<f>g, 

jA0).  For  i = 0,  1,...,N  and  j = 0,  1,...,M.  Also  , and  will 

denote  the  values  of  the  matrices  B,  C,  and  D,  respectively,  at  the  point 
0*>o+iA<J>,  jA0) . 

Equation  (35)  is  hyperbolic  if  there  is  a nonzero  crossflow.  The 
three  characteristics  at  a point  (<j>,9)  lie  between  the  angles  a and  a + 0; 
the  equation  is  parabolic  at  a point  if  0 is  zero,  as  the  characteristics 
have  the  same  tangent  or,  equivalently,  the  same  direction.  Along  a line 
of  flow  symmetry,  such  as  at  the  keel  or  at  the  waterline  on  a double 
model,  the  crossflow  is  zero  and  the  governing  Equation  (35)  is  parabolic. 
Consequently,  in  the  present  case  of  the  double  hull  models,  Equation  (35) 
is  a mixed  equation,  that  is  hyperbolic  and  parabolic  at  different  points 
of  the  region  of  integration  unless  the  crossflow  is  everywhere  zero.  In 
this  latter  case.  Equation  (35)  reduces  to  parabolic  form.  A solution 
method  which  is  applicable  to  both  parabolic  and  hyperbolic  equations  must, 
therefore,  be  used  to  solve  Equation  (35)  for  double  hull  models  or  hulls 

for  which  the  crossflow  is  zero  or  very  small  everywhere. 
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The  O'Brien  et  al. , implicit  finite  difference  scheme  is  used  to 
solve  Equation  (35).  It  is  a stable  scheme  for  any  positive  grid  spacing 
ratio  r = A0/A<f>  and  is  applicable  to  both  hyperbolic  and  parabolic 
equations.  It  consists  of  a one-step  forward  difference  in  the  ^-direction 
and  a central  difference  at  the  i + 1 step  in  the  9-direction.  In  this 
numerical  integration  scheme,  Equation  (35)  is  approximated  by  the  equation 


I 


D^i  ( W1+1J-W1J  \ B1 j /wi+1J+1-Wi+1^ -1 

hj,  V A<f>  / + h0  \ 2A0 


(38) 
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where  the  metrics  h*  and  h^  (defined  by  fc^“h^d$  and  fcg-hgdS)  are  evaluated 
at  the  point  ($q+1A$, JA0) . Rearranging  the  terms  of  Equation  (38) 
yields  the  equation 

- Wi+1  J"1  + DiJ  Wi+1J  + g«  W1+1  J+1 


C1J  - D1J  W1J  + h A<p  CiJ 
<P 


(39a) 


where 


’ft) 


(39b) 


The  values  of  W J for  j * 0,  1,...,M  are  specified  initial  conditions; 
they  may  be  obtained  from  experimental  data.  Along  the  symmetry  lines 
j ■ 0 (0-O)  and  j ■ M (0  ■-  y)  , the  crossflow  angle  8 is  zero,  so  that 
V*2  ■ t ■ 0,  the  equation  = 0 replaces  the  second  momentum  Equation  (36) 
along  the  two  lines  of  symmetry,  the  load  waterline  and  the  keel.  Moreover, 

91*22 

a - 0 on  these  lines,  so  that  g12>  g21»  g22»  h2»  h22’  S’  Cf  * A2*  and  3H 

ft 

are  identically  zero.  Thus,  along  the  symmetry  lines,  the  momentum  integral 
Equation  (3a)  reduces  to  the  equation 


39H  , n 3g12  9t  1 _ / H+2  9U 

3%  11  3t  9*0  2 % 11  V ♦ ' 0 9% 


and  the  auxiliary  Equation  (7)  reduces  to  the  equation 


0 nj1  + 9u  i + 6u  % ' r(H>  -6nG 


(1  a”  -*♦) 


(See  the  appendix  for  details).  The  system  of  Equation  (35)  can  be  used 
also  to  represent  Equations  (40)  and  (41)  together  with  t = 0,  if  the 
coefficient  matrices  B,  C,  and  D are  modified  to  the  following: 


and 


(42a) 


(42b) 


(42c) 


The  crossflow  angle  is  asymmetric  with  respect  to  0 - 0,  so  t is  also 
asymmetric  with  respect  to  this  line.  Accordingly,  the  central  difference 
approximation  for  3t/3f,Q  on  the  0 - 0 line  is 
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The  zero  crossflow  finite  difference  formula  analogous  to  Equation  (38)  Is 


After  rearranging  the  terms  of  Equation  (44),  one  obtains  the  equation 


where 


A similar  argument  at  the  symmetry  plane  0 »-  tt/2  yields  the  finite 
difference  equation 


(46a) 


The  finite  difference  Equations  (35),  (45),  and  (46)  form  a linear 
system  of  algebraic  equations  for  the  unknowns  Wi+1  ^ where  i - 0,1,... 
and  j - 0,1, . . . ,M. 

Let  matrix  D be  defined  by 


The  system  of  linear  equations  then  can  be  written  in  the  form 

DX  - C v (48a) 

where  X ■ (x^)C  and  C ■ (c^)  for  |i  ■ 1, . . . , 3(M+1),  where  the  components 
x^  of  the  vector  X are  defined  by 


X3k-Hc  " wk 


i+1  k 


(48b) 


where  K ■ 1,2,3  and  the  components  c of  the  vector  C are  defined  by 


c - Cik 
3k-Hc  tc 


(48c) 


Gaussian  reduction  is  used  to  solve  Equation  (48).  The  3x3  sub- 
matrices of  D have  been  inverted  explicitly  so  that  the  Gaussian  reduction 
of  Equation  (48)  is  very  fast  on  the  computer.  Back  substitution  is  used 
to  obtain  the  vector  X. 

COMPUTATIONAL  RESULTS  AND  DISCUSSION 
Some  sample  computational  results  of  the  boundary  layer  on  two  double 
ship  models  are  presented  in  this  section.  The  first  sample  computational 
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result  Is  the  boundary  layer  on  the  Lucy  Ashton  model  for  which  experl- 
mental  boundary  layer  data  are  given  by  Joubert  and  Matheson  and  which 
von  Kerczek*  computed  using  the  small  crossflow  approximation.  The  second 
sample  computational  result  described  below  is  the  boundary  layer  develop- 
ment on  the  Swedish  SSPA  Model  720  for  which  boundary  layer  experimental 
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data  and  calculations  are  provided  by  Larsson. 

It  is  first  necessary  to  describe  some  details  of  the  calculation  of 
the  surface  coordinate  system  before  embarking  on  a description  of  the 
results  of  the  boundary  layer  calculation.  There  are  many  different  surface 
coordinate  systems  that  one  can  use  for  three-dimensional  boundary  layer 
calculation  methods.  The  most  prevalent  coordinate  systems  used  for  ship 
boundary  layers  are  the  streamline  coordinate  system  and  the  coordinate 
system  made  up  of  the  cross-sectional  curves  and  their  orthogonal  tra- 
jectories on  the  hull  surface;  henceforth  referred  to  simply  as  the  cross- 
section  system.  The  streamline  coordinate  system*  has  the  advantage  of 
yielding  the  simplest  form  of  the  boundary  layer  equations,  but  it  may  be 
difficult  and  costly  to  generate  this  system  when  flows  about  a ship  hull 
at  nonzero  Froude  number  are  considered.  Thus  it  is  worthwhile  to  consider 
coordinate  systems  that  only  depend  on  the  ship  hull  geometry  and  not  on 
the  inviscid  flow. 

Figure  1 shows  a sample  of  the  cross-section  coordinate  system  recom- 
mended by  Mlloh  and  Patel^  on  the  Swedish  SSPA  Model  720.  The  calculation 
of  the  network  of  coordinate  lines  shown  in  Figure  1 is  described  in  the 
previous  section  of  this  report.  Calculation  of  the  Surface  Coordinate 
System  and  Potential  Flow.  The  main  feature  of  the  cross-sectional  co- 
ordinate system  that  has  been  found  to  be  objectionable  is  that  the  length- 
wise running  coordinate  lines  seem  to  diverge  on  certain  portions  of  the 
hull  (at  keel  near  the  bow  and  stern)  where  the  opposite,  i.e.,  convergence 
of  these  lines,  is  desirable.  Another,  minor,  annoyance  of  this  coordinate 
system  is  that  it  is  difficult  to  find  the  set  of  starting  values  at  any 
particular  station  for  the  coordinate  lines  along  the  length  (the  orthogonal 
trajectories  of  the  cross-sections)  that  will  result  in  a suitable  surface 
coordinate  grid.  Such  a grid  should  not  have  large  grid  intervals  or 


excessive  grid  crowding  on  some  portion  of  the  hull  (for  instance  note  the 

trajectories  on  the  bilge  and  and  those  at  the  load  waterline  near  the 

stern).  For  these  reasons  the  coordinate  system  shown  in  Figure  2 for 

SSPA  Model  720,  which  consists  of  the  lines  of  constant  0 and  their 

orthogonal  trajectories  of  constant  4>»  was  chosen  in  preference  to  the 

cross-section  system.  Recall  from  the  section  on  Calculation  of  the 

Surface  Coordinate  System  that  the  lines  of  constant  0 are  defined  by  the 

surface  representation  Equations  (23a, 23b).  Thus,  it  is  necessary  to 

always  first  represent  the  ship  hull  by  a surface  equation  of  the  type  of 

Equation  (23)  in  order  to  use  the  coordinate  system  of  Figure  2,  whereas 

the  cross-section  coordinate  system  does  not  require  a prior  analytical 

representation  of  the  ship  hull.*1  However,  the  calculation  of  the  ship 

surface  Equation  (23)  for  a typical  ship  hull  such  as  the  Swedish  SSPA 

Model  720  requires  only  about  one  to  one  and  a half  minutes  of  CDC  6700 

computer  execution  time.  This  calculation  of  the  surface  Equation  (23) 

(i.e.,  the  matrix  (A  ))  need  only  be  done  once  and  then  it  is  available  for 

mn 

several  other  uses.  Furthermore,  che  computer  method  used  to  calculate 

the  matrix  (A  ) is  an  old  one  and  several  modifications  of  this  method 
mn 

are  presently  under  development  that  are  expected  to  reduce  the  compu- 
tational time  by  a factor  of  about  100.  Thus,  the  need  to  calculate  the 
surface  representation  of  the  Equation  (23)  type  is  not  seen  as  a dis- 
advantage of  the  (0,$) -coordinate  system. 

The  Lucy  Ashton  double  model  boundary  layer  was  computed  using  the 
earlier  slender  body  theory  potential  flow  method^  because  the  first  test 
of  the  present  calculation  method  was  to  check  the  complete  crossflow 
formulation.  It  was  shown  previously  by  von  Kerczek  that  the  slender 
body  theory  gives  fairly  accurate  values  of  the  double-body  pressure 
distribution  on  the  Lucy  Ashton. 

The  boundary  layer  calculation  method  of  this  report  is  implemented 
in  terms  of  the  (0,4>)-coordinate  system  but  the  computed  boundary  layer 
results  are  given  in  terms  of  the  streamline  momentum  thickness  0^,  dis- 
placement thickness  6^,  shape  factor  H,  and  the  stream  coefficient  of  skin 
friction  C^.  This  is  done  to  facilitate  the  comparison  of  the  present 
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results  with  previous  boundary  layer  calculation  methods  and  experi- 
1 21 

ments.  ’ Earlier  calculations  of  the  Lucy  Ashton  double  body  boundary 
layer  by  von  Kerczek  showed  that  crossflow  is  small  almost  everywhere  on 
the  hull. 

Figure  3 shows  several  streamlines  computed  by  slender  body  potential 
flow  theory.  In  bottom,  elevation,  front,  and  rear  views  on  the  Lucy 
Ashton  double  model.  Computed  boundary  layer  results  will  be  shown  along 
these  streamlines.  Figures  4a-4d  show  a comparison  of  the  distribution  of 
the  streamline  momentum  thickness  9^,  along  the  streamlines  1,  6,  10,  and 
13  shown  in  Figure  3b,  as  computed  by  the  present  complete  crossflow 
method  and  the  small  crossflow  method  of  von  Kerczek.*-  Figures  5a-5d  show 
the  distribution  of  streamline  skin  friction  coefficient  and  Figures  ba- 
be show  the  distribution  of  the  crossflow  angle  in  radians  along  these  same 
streamlines.  Note  from  Figures  3 through  b that  there  is  little  difference 
in  the  boundary  layer  characteristics  that  are  predicted  by  the  present 
complete  crossflow  method  and  the  small  crossflow  method.  This  is  not  an 
unexpected  result  because  the  Lucy  Ashton  is  a fairly  slender  hull  with 
very  slowly  changing  cross-section  shape  along  the  length  of  the  ship. 

Hence  the  values  of  the  coefficient  ^ are  small  everywhere  along  the  hull 
and  it  is  reasonable  to  expect  fairly  small  boundary  layer  crossflow 
effects.  The  differences  in  the  two  sets  of  results  shown  in  Figures  3 
through  b are  due  mainly  to  the  differences  in  the  numerical  integration 
method  used  by  von  Kerczek-*  and  the  present  method.  This  is  indicated  by 
the  differences  in  the  results  on  the  keel,  shown  in  Figures  4d  and  5d, 
where  the  two  methods  solve  identical  equations. 

The  second  test  calculation  is  of  the  boundary  layer  on  the  Swedish 

21 

SSPA  Model  720  double  body.  Figure  7,  taken  from  Larsson’s  report,  shows 
front  and  rear  views  of  the  streamlines  along  which  measured  and  computed 
boundary  layer  properties  were  given.  Larsson’s  calculation  method  starts 
from  a momentum-integral-entrainment  method  closely  related  to  the  one 
described  in  this  report.  The  main  difference  between  these  two  boundary 
layer  calculation  methods  is  in  the  auxiliary  data  used  for  the  crossflow 

velocity  profiles  and  the  numerical  implementation  of  the  methods. 

21 

Larsson's  method  computes  the  boundary  layer  in  the  streamline  surface 
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Bottom  and  Side  Views  of  Streamlines 


Figure  3 (Continued) 
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- Streamline  Momentum  Thickness  versus  Axial  Distance  for  Flow 
Along  Streamlines  of  the  Lucy  Ashton 
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Figure  4 (Continued) 
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Figure  4b  - Streamline  6 


Figure  5 (Continued) 
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Figure  6 - Crossflow  Angle  versus  Axial  Distance  for  Flow  Along 
Streamlines  of  the  Lucy  Ashton 
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Figure  6a  - Streamline  1 
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Figure  6b  - Streamline  6 
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Figure  6c  - Streamline  10 


coordinate  system  and  uses  a numerical  integration  method  that  includes 
the  crossflow  terms  by  an  explicit  method.  The  present  calculation 
scheme  may  improve  on  Larsson's  results  only  in  detail  but  not  in  a funda- 
mental way. 

The  surface  coordinate  grid  on  which  the  present  calculations  were 
made  is  similar  to  but  slightly  coarser  than  the  one  shown  in  Figure  2. 

The  grid  spacing  is  A0  * 6 deg  and  As  ■ 0.025  (As  ^ A<J>).  The  Chang 
potential  flow  method  was  used  to  compute  the  values  of  the  potential  at 
the  center  of  rectangular  elements  in  the  (6,$)  plane.  These  potential 
flow  rectangular  elements  (A0-1O  deg,  A4>«0.045)  were  much  larger  than  the 
grid  elements  of  the  boundary  layer  calculation  so  that  the  interpolation 
technique  with  the  cubic  spline-on-spline  described  in  the  section  on  The 
Ship  Surface  Coordinate  System  was  used. 

Initial  conditions  for  the  quantities  0^,  H,  and  t - tan  0 at  station 

s “ -0.5  (which  very  nearly  coincides  with  one  of  the  ©-coordinate  curves 

21 

on  the  surface  of  the  hull)  were  obtained  from  Larsson's  experimental 

results.  These  Initial  data  are  shown  in  Figure  8.  Intermediate  values  of 

the  data  shown  in  Figure  8 were  obtained  by  linear  interpolation  because 

it  was  felt  that  the  sparsity  and  quality  of  the  experimental  data  did  not 

justify  a more  accurate  interpolation.  Figures  9 through  13  show  the 

distributions  of  the  streamline  momentum  thickness  0^,  crossflow  angle  3, 

and  skin  friction  coefficient  C^  along  the  streamlines  1 through  8 of 

Figure  7.  In  each  case,  the  results  of  the  present  calculation  are 

21 

compared  to  the  corresponding  results  of  Larsson's  experiment  and  his 
calculation  that  includes  the  complete  crossflow  but  not  the  modification 
of  the  hull  offsets  by  the  values  of  the  local  displacement  thickness.  It 
can  be  seen  in  Figures  9 through  13  that  the  present  predictions  seem  to 
correlate  with  the  experimental  data,  on  the  average,  about  as  well,  or 
possibly  slightly  better  than,  Larsson's  computational  results.  In 
particular,  the  present  boundary  layer  predictions  on  streamline  5 extend 
nearly  to  the  stern  of  the  model,  in  fairly  good  agreement  with  the  experi- 
mental data,  whereas  Larsson's  results  on  this  streamline  seem  to  terminate 
somewhat  earlier,  apparently  because  of  some  breakdown  in  his  calculation 
method. 


U2 


Figure  7 - Streamlines  on  the  SSPA-720  Model  Along  Which  Larsson  Measured 

the  Boundary  Layer 


Figure  8 - Starting  Values  of  H,  $ and  0..  at  Station  S ■-  0.5  on 

SSPA-720  Model 
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Figure  9 - Boundary  Layer  Characteristics  versus  Axial  Distance 
for  Flow  Along  Streamline  1 on  Model  SSPA-720 

X 103 


s 


Figure  9a  - Streamline  Momentum  Thickness 
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Figure  9 (Continued) 
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Figure  9b  - Crossflow  Angle  8 


Figure  10  - Boundary  Layer  Characteristics  versus  Axial  Distance  for  Flow 
Along  Streamline  3 on  Model  SSPA-720 


Figure  10  (Continued) 
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Figure  10b  - Crossflow  Angle  B 
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Figure  11  (Continued) 


Figure  lib  - Crossflow  Angle  B 
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Figure  12  (Continued) 


Figure  12b  - Crossflow  Angle  8 
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Figure  13  - Boundary  Layer  Characteristics  versus  Axial  Distance  for  Flow 
Along  Streamline  8 on  Model  SSPA-720 
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gure  13  (Continued) 
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Figure  13b  - Crossflow  Angle  8 


It  is  to  be  noted  that  the  experimental  values  of  the  crossflow  angle 
3 shown  in  Figures  9b,  10b,  lib,  12b,  and  13b  all  exhibit  a change  in  sign 
(hence  crossflow  reversal)  between  the  stations  s * 0.6  and  s * 0,8.  The 
results  of  the  boundary  layer  calculations  shown  in  these  figures  consist- 
ently miss  this  flow  reversal.  This  indicates  that  the  crossflow  model  may 
need  considerable  improvement  in  order  to  reliably  predict  crossflow 
reversal  as  the  boundary  layer  approaches  separation.  However,  the  cross- 
flow  discrepancy  at  the  stern  of  the  SSPM  Model  720  may  also  be  due,  in 
large  part,  to  the  discrepancy  between  the  potential  flow  pressure 
distribution  and  the  actual  pressure  distribution  at  the  stern.  The 
experimental  and  computational  results  indicate  that  the  boundary  layer 
is  very  thick  at  the  stern  of  the  model;  hence,  the  pressure  distribution 
must  be  different  from  the  potential  flow  values  there.  Note  also  that  the 
degree  of  accuracy  in  predicting  the  primary  quantities  of  interest,  the 
boundary  layer  momentum  thickness  0^  and  skin  friction  coefficient  C^,  is 
considerably  better  than  the  prediction  of  the  relatively  small  values  of 
the  crossflow  angle  0. 

On  the  basis  of  the  comparisons  with  experimental  boundary  layer  data 
shown  in  Figures  9 through  13,  the  overall  assessment  of  the  present 
boundary  layer  calculation  method  is  that  it  can  predict  boundary  layers  on 
relatively  fine  double  ship  models  with  fair  accuracy  to  within  a distance 
of  the  stern  of  about  10  percent  of  the  ship's  length.  In  this  area,  the 
boundary  layer  thickens  very  rapidly  and  approaches  separation.  Calcu- 
lation of  this  near-separated  boundary  layer  region  must  await  further 
developments  of  boundary  layer  theory. 

CONCLUDING  REMARKS 

This  report  presents  a momentum  integral  method  for  computing  three- 
dimensional  boundary  layers  for  ships.  Most  of  the  technical  details  for 
carrying  out  the  computational  problem  of  solving  the  momentum- integral 
boundary  layer  equations  are  worked  out  here  and  have  been  implemented  in  a 
set  of  computer  programs.  The  basic  method  can  be  used  to  calculate 


certain  boundary  layer  parameters,  such  as  boundary  layer  thickness  or 
skin  friction,  with  fair  accuracy  over  a large  portion  of  hulls  where 
unseparated  flow  is  maintained.  The  computer  programs  are  now  ready  to  be 
modified  so  as  to  improve  the  crossflow  modeling,  with  whatever  new 
experimental  data  becomes  available.  Alternatively,  portions  of  the 
developments  described  in  this  report,  such  as  the  surface  coordinate 
system  and  inviscid  flow  calculation,  can  be  used  in  other  methods  for  cal- 
culating the  boundary  layer. 
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APPENDIX 

MATHEMATICAL  DETAILS 


Substituting  directly  from  Equation  (21)  into  Equation  (6)  yields 
011  = 911  (cos2a  _ t(f1+f2)  sin  a cos  a + t2f3  sin2a]  = 0llgll(t,H)  (49a) 
0^2  = ®n  sin  a cos  a + cos2a  “ tf ^ sin2a]  = O^g^Ct.H)  (49b) 

®2i  = ®n  sin  a cos  a + tf^  cos2a  - tf2  sin2a]  = ®n821^t,H^  (49c) 

022  = 011  fsin2a  + sin  a cos  a + t2f3  cos2a]  = 0ug22(t,H)  (49d) 


A1 

= 0, , (H  cos  a - tf , 
11  4 

sin  a)  = O^h^t.H) 

(49e) 

A2 

= 0, , (H  sin  a + tf . 
11  4 

cos  a)  = 01;Lh2(t,H) 

(49f) 

The  momentum  integral  Equations  (3)  can  then  be  written 
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where 
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Use  is  made  here  of  the  identity 
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Let  T and  S be  defined  by 
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The  entrainment  Equation  (7)  can  then  be  written  in  the  form: 
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Moreover,  using  Equations  (10),  (49e),  and  (49f) 
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The  final  form  of  the  entrainment  equation  is 
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ONR/Pasadena 

1 

NAVSHIPYD  CHASN/Lib 

ONR/San  Francisco 

1 

NAVSHIPYD  LBEACH/Lib 

NORDA 

2 

NAVSHIPYD  MARE 

USNA 

1 Library 

1 Code  250 

1 Tech  Lib 

1 Nav  Sys  Eng  Dept 

1 

NAVSHIPYD  PUGET/Lib 

1 B.  Johnson 

1 

NAVSHIPYD  PEARL/Code  202.32 

NAVPGSCOL 

1 Library 

8 

NAVSEC 

1 T.  Sarpkaya 

1 SEC  6034B 

1 J.  Miller 

1 SEC  6110 

NADC 

1 SEC  6114H 

1 SEC  6120 

NOSC 

1 SEC  6136 

1 SEC  6140B 

1 Library 

1 SEC  6144 

1 T . Lang 

1 SEC  6148 

1 J.W.  Hoyt 

1 D.M.  Nelson 

1 

NAVSEC,  NORVA/ 6660 . 03 , Blount 

NCSC/712,  D.  Humphreys 

12 

DDC 

NCEL/ Code  131 

1 

AFOSR/NAM 

NSWC,  White  Oak/Lib 

1 

AFFOL/FYS,  J.  Olsen 

Copies 


Copies 


2 


1 

3 


1 

1 

1 

2 


1 

4 


2 

5 


1 

1 

1 


MARAD 

1 Div  of  Ship  R&D 
1 Lib 

NASA  HQS/Lib 
NBS 

1 Lib 

1 P.S.  Klebanoff 
1 G.  Kulin 

NSF/Eng  Lib 

LC/Sci  & Tech 

DOT/Lib  TAD-491.1 

MMA 

1 National  Maritime 
Research  Center 
1 Library 

U.  of  Bridgeport/E.  Uram 

U.  of  Cal/Dept  Naval  Arch, 
Berkeley 
1 Library 
1 W.  Webster 
1 J.  Paul ling 
1 J . Wehausen 

U.  of  Cal,  San  Diego 
1 A.T.  Ellis 
1 Scripps  Inst  Lib 

CIT 

1 Aero  Lib 
1 T.Y.  Wu 
1 A. J.  Acosta 
1 I.  Saber sky 
1 D.  Coles 

City  College,  Wave  Hill/ 
Pierson 

Catholic  U.  of  Amer. /Civil  & 
Mech  Eng 

Colorado  State  U./Eng  Res  Cen 


1 U.  of  Connecticut/Scottron 

1 Cornell  U./Shen 

2 Florida  Atlantic  U. 

1 Tech  Lib 
1 S . Dunne 

2 Harvard  U . 

1 G.  Carrier 
1 Gordon  McKay  Lib 

1 U.  of  Hawaii/Bretschneider 

1 U.  of  Illinois/ J.  Robertson 

4 U.  of  Iowa 

1 Library 
1 L . Landweber 
1 J . Kennedy 
1 V.C.  Patel 

1 Johns  Hopkins  U. /Phillips 

1 Kansas  State  U. /Nesmith 

1 U.  of  Kansas /Civil  Eng  Lib 
1 Lehigh  U. /Fritz  Eng  Lab  Lib 

5 MIT 

1 Library 
1 P.  Leehey 
1 P.  Mandel 
1 M.  Abkowitz 
1 J . N . Newman 

4 U.  of  Minn/St.  Anthony  Falls 

1 Silberman 
1 Lib 
1 Song 
1 R.  Arndt 

4 U.  of  Mich/NAME 

1 Library 
1 F.  Ogilivie 
1 Hammitt 
1 Cough 
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VP  I 

1 H.L.  Moses,  Dept  Mech 
Eng 

1 D.P.  Telionis,  Dept  Mech 
Eng 

1 J.  Schetz,  Dept  Aero  & 
Ocean  Eng 

Webb  Inst 
1 Library 
1 Lewis 
1 Ward 

Woods  Hole/Ocean  Eng 
Worchester  PI/Tech  Lib 
SNAME/Tech  Lib 

Bethlehem  Steel/Sparrows  Point 

Bethlehem  Steel/New  York/Lib 

Bolt,  Beranek  & Newman/Lib 

Exxon,  NY/Design  Div,  Tank 
Dept 

Exxon  Math  & System,  Inc. 

General  Dynamics,  EB/ 
Boatwright 

Gibbs  & Cox/Tech  Info 

Hydronautics 
1 Library 
1 E.  Miller 
1 V.  Johnson 
1 C.C.  Hsu 

Lockheed,  Sunnyvale/Waid 

McDonnell  Douglas,  Long  Beach 
1 T.  Cebeci 
1 J.  Hess 

Newport  News  Shipbuilding/Lib 
Nielsen  Eng  & Res 


Copies 

Copies 

Code 

Name 

1 

Oceanics 

1 

1544 

R.  Cumming 

1 

1544 

R.  Boswell 

3 

Rand 

Corp 

1 

1544 

E.  Caster 

1 

E.R.  Van  Driest 

1 

1544 

S.  Jessup 

1 

C.  Gazley 

1 

1552 

T . T . Huang 

1 

J.  Aroesty 

1 

1552 

N.  Salvesen 

30 

1552 

C.  von  Kerczek 

1 

Rockwell  International/B. 

1 

1552 

H.T.  Wang 

U j ihara 

1 

1552 

J.  McCarthy 

1 

1552 

T.  Langan 

1 

Sperry  Rand/Tech  Lib 

1 

1560 

G.  Hagen 

1 

Sun  Shipbuilding/Chief  Naval 

1 

1562 

M.  Martin 

Arch 

1 

1564 

J.  Feldman 

1 

Robert  Taggart 

1 

1568 

G.  Cox 

1 

Tracor 

1 

1572 

M.D.  Ochi 

1 

1572 

E.  Zarnick 

2 

Westinghouse  Electric 

1 

1572 

C.M.  Lee 

1 

M.S.  Macovsky 

1 

1576 

W.E.  Smith 

1 

Gulino 

1 

1615 

R.J.  Furey 

CENTER  DISTRIBUTION 

1 

184 

H.  Lugt 

Copies 

Code 

Name 

1 

1843 

J.  Schot 

1 

1843 

C . Dawson 

1 

11 

W.M.  Ellsworth 

1 

117 

R.M.  Stevens 

1 

19 

M.M.  Sevik 

1 

1942 

J.T.  Shen 

1 

1500 

1 

1942 

W.R.  Brown 

1 

1504 

V.J.  Monacella 

1 

1942 

F.C.  DeMetz 

1 

1942 

T.M.  Farabee 

1 

1506 

M.K.  Ochi 

1 

1942 

F.E.  Geib 

1 

1507 

D . Moran 

1 

1942 

T.C.  Mathews 

1 

1508 

F.  Peterson 

10 

5211.1 

Reports  Distribution 

1 

1521 

P.  Pien 

1 

522.1 

Unclassified 

Lib 

(C) 

1 

1524 

Y.T.  Shen 

1 

522.2 

Unclassified 

Lib 

(A) 

1 

1524 

W.C.  Lin 

1 

1524 

W.  Day 

1 

1532 

G.  Dobay 

1 

1532 

M.  Wilson 

1 

1540 

W.B.  Morgan 

1 

1541 

P.  Granville 

1 

1542 

B.  Yim 

